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Specialization of Appell's functions to univariate 
hypergeometric functions 



Raimundas Vidunas 
Kobe University* 



Abstract 

Univariate specializations of Appell's hypergeometric functions Fi, F2, -F3, F4 
satisfy ordinary Fuchsian equations of order at most 4. In special cases, these dif- 
ferential equations are of order 2 and could be simple (puUback) transformations 
of Euler's differential equation for the Gauss hypergeometric function. The pa- 
per classifies these cases, and presents corresponding relations between univariate 
specializations of Appell's functions and univariate hypergeometric functions. The 
computational aspect and interesting identities are discussed. 

1 Introduction 

Appell's hypergeometric functions Fi{x, y), ^2(2;, y), F'i{x, y), Fi^x, y) are defined by the 
following double hypergeometric series: 



Fi 



a; bi,b2 
c 

a; 61,62 



Ci , C2 

01,02; 61,62 



Fa 



a; 6 

Cl,C2 



X, y 



X, y 



X, y 



X, y 



00 00 



EV^ {a)n+m (61 )n (62)ni „ ra 



71^0 771 — 

00 00 



{c)n+,n n\m\ 



Y^ V^ {a)n+m (6l)n (62)m „ „ 

00 00 

EE 

n— m—0 

00 00 



(ai)n(a2)m(6l)„(62)m „ ,„ 

X y , 



{c)n+m n\ m 



{a)n+m (6) 



n n (Cl)„ (C2)m n\ m\ 



n+m ^nyra_ 



They are bivariate generalizations of the Gauss hypergeometric series 

A, B 



2F1 



C 



^ {A)n{B)n r 



(1) 

(2) 
(3) 
(4) 

(5) 
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In particular, Appell's functions satisfy Fuchsian (that is, linear homogeneous and with 
regular singularities) systems of partial differential equations, holonomic of rank 3 or 4, 
which are analogous to Euler's hypergeometric equation for the 2F1 function. Euler's 
equation for ([5]) is 

z(l - z) ^-0- + {C-{A + B + l)z) ^ -AB y{z) = 0. (6) 

It is a Fuchsian equation with three singularities: z — 0, z — I and z = cx). Local 
exponent differences are equal, respectively, to 1~C,C — A — B,A~B a.t them. 

When the two arguments x, y of Appell's functions are algebraically related, the 
univariate specializations satisfy Fuchsian ordinary differential equations of order at most 
4. This paper considers the following questions: 

Which univariate specializations of Appell's functions satisfy a second order 
Fuchsian ordinary differential equation? 

Which univariate specializations of Appell's functions satisfy Euler's equa- 
tion ^ for some parameters A, B, C, up to a projective transformation 
y{z) <—> 9{z)y{z) or a more general pull-back transformation 

z^ip{x), y{z)v-^e{x)y{ip{x)), (7) 

where 6{z) or 9{x) is a power function, and 1^9(2;) is a rational function? 

Other interesting questions are: Which univariate specializations of Appell's functions 
satisfy a Fuchsian ordinary differential equation of order 3? Up to pull-back transforma- 
tions, which univariate specializations of Appell's functions satisfy differential equations 
for generalized hypergeometric functions 3F2(z) or 4F3(z)? 

The paper classifies univariate specializations of Appell's -Pj, -F3, -F4 functions satis- 
fying second order ordinary Fuchsian equations, when those ordinary equations follow 
from the partial differential equations for Appell's functions without reductions due to 
factorization of the respective differential operators (or reducibility of the monodromy 
group). In particular, we do not consider isolated solutions of reducible Appell's systems 
of partial differential equations that are expressible via Gauss hypergeometric functions. 
Definition 11.11 below specifies the way ordinary differential equations follow fully from 
systems of partial differential equations. 

We identify those univariate specialization of Appell's F2, F-^, F4 functions that can 
be expressed in terms of Gauss' hypergeometric function ^. Their ordinary differential 
equations are pull-back transformations of Euler's equation ([S]), and follow fully from 
the respective partial differential equations. To deal with univariate specializations onto 
singularity curves of the respective systems of partial differential equations, we present 
general ordinary Fuchsian equations for them. In the cases of Fi, F2, F3 functions, the 
specializations onto singularity curves can be expressed in terms of 2F1 or 3F2 hypergeo- 
metric functions. 

In Appendix, we briefly review identities between Appell's and Gauss hypergeometric 
functions existing in literature. Our explicit results that appear more or less novel are 
the following: 



• A case of F2{x,2 — x) function that can be expressed as a 2F1 function with a 
quadratic argument; see Theorem l2.4l Equivalent identification of ^3(2;, x/{2x — 1)) 
functions was noticed by Karlsson 9 . These relations apply nicely when the nion- 
odromy group of the 2F1 function is dihedral; then the F2 and F3 solutions of the 
same ordinary Fuchsian equation are terminating, and dihedral 2F1 functions can 
be expressed as elementary functions. We show this application in Section [71 

• A separation of variables case for the F2 function; see formula ([55)1 . It is related to 
the well-known Bailey case [4] of variable separation for the F4 function 



F4 



c, a+6 — c+1 



xil-y),y{l -x) 



2F1 



a, b 
c 



2F1 



a, b 
a+6-c+l 



y] (8) 



via a known transformation between F2 and F/^ functions; 

• Translation of the above two F2 cases to relations of Fi and F3 functions to Gauss 
hypergeometric functions. 

• A few cases of -^^(i^, (1 — t)'^) functions (that is, specializations to the quadratic 
singular curve of F4) expressible as 2F1 or 3F2 functions, or products of two 2F1 
solutions of the same Euler's equation © . 

Identities between bivariate and univariate hypergeometric series are usually derived 
by methods of series manipulation or using integral representations. The method of re- 
lating hypergeometric functions as solutions of coinciding differential equations is usually 
considered tedious and computationally costly. With powerful computer algebra tech- 
niques available, the method of identifying differential equations can be worked out rather 
comprehensively. We discuss general computational techniques in Section [S] Here below 
and in Section [5] we gradually introduce the algebraic setting, computational methods 
and some shortcuts. 

Consider a holomorphic function F{x,y) on an open set in C x C, and a univariate 
specialization F{x{t),y{t)) of it. The full derivatives with respect to t are expressed 
linearly in terms of the partial derivatives dF/dx, dF/dy, d^F/dx^, etc. For example. 



dF 
Itt 
(fF 



dF 
dx 

dx"^ 



y 



dF 
dy' 

2xy 



d^F 
dxdy 



+r 



d^F 
dy^ 



dF .. dF 



dx 



y 



dy' 



(9) 
(10) 



Here x = dx/dt, y — dy/dt, x — d^x/dt^, y = d^y/dt^. To compute next order full 
derivatives, one applies the Leibniz rule and lets d/dt act on the partial derivatives by 
copying the action on F in ^. 

We identify partial differential equations with ordinary differential equations in the 
specialization variable t by identifying their "mixed" forms in the partial derivatives 
with the coefficients specialized to functions in t. To fix a univariate specialization 
mapping, consider a holomorphic map <& from an open subset of C to C x C, mapping 



t 1-^ ix{t), y{t)). If F{x, y) is a holomorphic function on the image of $, then F o <J) is 
the univariate speciahzation F(x{t),y{t)). Geometrically, a univariate specialization is 
the puUback of F{x, y) with respect to $. 

Defintion 1.1 If _E is an ordinary differential equation with respect to t, its partial 
differential form under $ is the expression where the derivatives with respect to t are 
replaced by partial derivatives following formulas (|9|)- (fT0l) . etc. 

If .E is a partial differential equation with respect to x, y, its specialized form under $ 
is the expression where the coefficients to the partial derivatives are specialized x i— *■ a;(t), 

y^ y{t)- 

Let H denote a system of partial differential equations with respect to x, y. The 
ordinary differential equation E is said to follow fully from H under $, if its partial 
differential form under <I> coincides with the specialized form (under $) of some partial 
differential equation following from H by algebraic and partial differentiation operations. 

We also say that E is implied fully by H, if there is a specialization map $ such that 
E follows fully from H under $. Here we allow H, E and $ to have parameters; the 
parameters of H may specialize (to constants or functions in the parameters of E, $) in 
the specialized forms of its partial differential equations under <I>. 

The classical systems of partial differential equations for AppeU's functions F2, F3, 
F4, Fi are presented, respectively, in formulas (HU-IHD, §^~^, 1^-^, (|7n)) -(|7T |) 
below. Algebraically, linear partial differential equations with polynomial coefficients are 
identified with partial differential operators in the Weyl algebra C[x,y]{d/dx,d/dy). A 
system of partial differential equations corresponds to a left ideal in the Weyl algebra. 

The systems of partial differential equations for F2, F3, F4 functions have rank 4. 
Correspondingly, Grobner bases for the mentioned left ideals in the Weyl algebra give 
linear expressions of all partial derivatives in terms of 4 partial derivatives of low order, 
say d^F/dx^ dF/dx, dF/dy, F. If we take a total degree ordering of partial derivatives, 
the leading coefhcients in the expressions for higher order derivatives vanish only when 
specialized onto a singularity curve of the differential system. 

If a univariate specialization of AppeU's F2, F3 or F^ function is not onto a singularity 
curve of the differential system, the full derivatives can be expressed linearly in terms of 
the 4 basic partial derivatives, and a linear relation between d'^F/dt'^, d'^F/dt^, d^F/dt^ , 
dF/dt, F gives an ordinary differential equation of order at most 4 for the univariate 
function. If the univariate specialization is onto a singularity curve, it turns out that we 
get an ordinary differential equation of order less than 4. All these ordinary differential 
equations are Fuchsian. 

Similarly, the system of partial differential equations for the Fi function has rank 3, 
and its univariate specializations satisfy ordinary Fuchsian equations of order at most 3. 

We are interested in cases when the ordinary differential equation has order 2. Such 
a second order differential equation must follow from partial differential equations of 
order at most 2. For non-singular specializations oi F2, F3, F4 functions, linear relations 
between partial derivatives of order at most 2 are generated by the classical pairs of 
partial differential equations. 



2 Identities with Appell's F2 function 



A system partial differential equations for the F2{x, y) function is: 

Q2p Q2p gp Qp 

x(l - x)—-^ - xy——- + (ci - (a + 61 + l)x) biy- abiF = 0, (11) 

,it,^ ijxoy ox ay 



dx"^ 



dx 

u ±' d'^F dF dF 

y(l - y^in - ^yinr + (^2 -ia + b2 + l)y) ^- - b^x^- - ab2F = 0. (12) 

oy^ oxoy oy ox 

This is an equation system of rank 4. Its singular locus on P^ x P^ is the union of the 
following lines: 



X = Q, X — 1, X — oo, J/ = 0, J/ = 1, y — oo, X + y = \. 



(13) 



First we handle univariate specializations onto the singular lines. We determine that 
generally these specializations satisfy hypergeometric equations for 2F1 or 3F2 functions. 
Explicit identification of the following three function pairs is presented after Theorem 
Those identities can be easily derived from power series manipulations as well. 



Theorem 2.1 The following function pairs satisfy the same ordinary differential equa- 
tions (of order 2 or 3/' 



• Fr 



• Fr 



• Fo 



a; 61,62 

Cl,C2 

a; 61,62 

Cl,C2 

a; 61,62 

Cl,C2 



X, I and 2F1 



a, 61 

Cl 



a, 61, a — C2 + 1 



X, 1 j and 3F2 

x,l — X \ and (1 — x)^" 3F2 



Cl, a + 62 - C2 + 1 

a, Cl - 61, a- C2 + 1 
Cl , a + 62 - C2 + 1 



X- 1 



Proof. The first specialization is trivial: after the substitution y = into the coefficients 
of pip we have differentiation with respect to x only. 

To handle the other two cases, we use the third order differential equation for the 
3F2(a:) function; see formula pOOp in the Appendix. 

For the ^2(2:, 1) function, subtract (I12p from pT|) and set y = 1 in the coefhcients: 



d'^F dF 

x{l - x)-^ + (ci - (a + 61 - 62 + 1)2:) — 

dF 
-(c2-a + 6i-62-l)^— -a(6i-62)F = 0. 
oy 



(14) 



Then make the following combination of partial differential equations for the F2{x, 1): 



a;— [Eq. ([T4D] + (a + 62-C2 + l)[Eq. ([TTD] - 61 [Eq. IHD]. 



(15) 



After setting y = 1 in the coefficients, we get an equation with differentiation by x 
only. That differential equation is identified as hypergeometric equation (jlOOp with the 
indicated values of the parameters A, B, C, D, E. 



For the F2{x, 1 — x) function, add pT|) and P^ and set y = 1 ~ x: 

/Q2p Q2p g2p\ gp 

dF 

+ {c2~a-bi-b2-l + ia + bi+b2 + l)x)^--a{bi+b2)F = 0, (16) 

oy 

and consider the following combination of partial differential equations for it: 

jL _ I , ^) p,. „ , ^^^^[E,. ODl . ^lE. 031. ,17, 

With y — 1 — X, we recognize the full derivatives 

d^F _ d^F d^F d^F d^F 

dx^ dx^ dx'^dy dxdy"^ dy^ ' 

£p _ d^F d^F d^F dF _ dF dF 

dx^ dx"^ dxdy dy^ ' dx dx dy 

in ()17|) . and get an ordinary differential equation of order 3, with the following singular- 
ities and local exponents: 

at a; = : 0, 1 — ci, C2 — a — 62; 
at a; = 1 : 0, 1 — C2, ci — a — fei; 
at X = 00 : a, a+l, 61+62- 

The differential equation for the 3F2 function (times the power factor) has the same data, 
making the requisite check of the differential equation on the second function worthwhile. 

□ 

The cases when 3F2 functions in Theorem 12.11 become 2F1 functions can be found by 
equating a pair of upper and lower parameters. Note that the case C2 =0 — 61 + 62 + 1 for 
the i^2(a;, 1) function is already evident from (|14p . while the case ci + C2 = a + 61 + 62 + 1 
for the F2{x, 1 — x) function is visible in (|16p . Euler's equation ([6]) is easily recognizable 
in those two formulas. 

The next two theorems tell that second order Fuchsian equations for univariate F2 spe- 
cializations outside the singularity lines, implied fully by (fTT1) - P^ . follow from pT|) - p^ 
without differentiation. 

Theorem 2.2 The following univariate specializations of Appell's F2{x,y) function sat- 
isfy ordinary Fuchsian equations of second order (with respect to x or t): 



'' 261,262 



X, 2 - X , F2 



bi+b2- h 61,6 



■J2 



261,262 
In the second function, s is any constant. 



.2 , (t 



1-^'1-V^ ■ (20) 



Proof. For the F2{x,2 — x) function, divide pT|) by x and (|12p by —y, add the two 
equations, and substitute y = 2 — x in the coefficients: 

, fd^F „ d^F d^F\ fci , b2X \ dF 



\ dx'^ dxdy dy^ J \x x ~ 2 J dx 

"' • a + b2 + l-b,^-^\'-l-a('-l + J^\F = 0. (21) 



x— 2 X J dy \x x—2 

If ci = 26i and C2 = 262, foUowing (fT9| we recognize the second order equation 

,d2^ /26i 262 \ dF fbi 62 \ 

ax \ X x—2 J ax \x x — 2 / 

For the second F2 function, the foUowing differential equation can be checked by sub- 
stituting (P))- (jlOp with evaluated i, y, i, j/, then linearly eliminating the derivatives 
d'^Fjdx^, d^Fjdy^ using equations pT|) ~ (fT^ : 

d^F f Abit 452(t + s) \dF , ^. f 46i Ah 



iH^-('^-'-«7^-FtSttK^°^'-) 



dt^ \t^~l t^ + 2st+lj dt ^ "^'yt^-l t^ + 2st+l^ 

a 

Theorem 2.3 Suppose that a univariate specialization of Appell's F2{x,y) function sat- 
isfies a second order ordinary Fuchsian equation fully implied by the system (|ll|) - (|12p . 
Then either the specialization is into singular locus (|13p . or they are represented by one 
of the two functions in (|20p . 

Proof. We can assume that x{x ~ \)y{y ~ l){x + y ~ \) ^ Q. If F{t) = F2 {x{t),y{t)) 
satisfies a second order Fuchsian equation, second order partial derivatives come from 
the d'^F/dt^ term following ([10]). The rank of the differential system (ITT])-(Il2l) is 4, and 
there are 2 linearly independent relations between partial derivatives of order at most 2. 
Therefore equations pT|) - P^ linearly generate all partial differential equations of order 
2 that algebraically follow from them. (In other words, the corresponding two differential 
operators in the Weyl algebra C[x, y] <d/dx, d/dy> linearly generate order 2 operators 
of the corresponding left ideal.) 

Linear elimination of d^F/dx^, d^F/dy^ from (jlOp using equations (fTT1) - P^ gives 
the following coefficient to d^F/dxdy. 

-J^x^ + 2xy+-^f. (24) 

I — X 1 — y 

This expression must be equated to 0. The expression is a quadratic form in the deriva- 
tives x,y, with the discriminant equal to 4(1 — x — y) / (l — x){l — y) . We can factorize (|24|) 
into linear differential forms if we parametrize the surface {1 — x — y)u^ — (1 — a;)(l — y). 
Accordingly, we substitute 

.^ "-5-/' ^ (25) 



The parameter u is undefined only ii x + y = 1. 
Expression ((24)) then factorizes as follows: 



(2x{x — l)du — {u — l){x — u — l)dx)(2x{x ~ l)du — {u + l){x + u — l)dx) 

(x + u^-lf ' 

where we write du, dx instead of du/dt, dx/dt, because the variable t is irrelevant. (It can 
be composed with any univariate function without changing our differential relations.) 
The difference in the two numerator factors is in the sign of u. The differential equation 

2x{x - l)du ^ (u - l){x - u - l)dx (27) 

is generally solved by the following algebraic relation between u and x: 

{u + x-lf + Cix-l){u-lf = 0, (28) 



where C is an integration constant. After eliminating u from (p5)l and (|28p we get the 
algebraic relation between x and y that solves ([M)) : 



(x + y - 2)2 - 2C{y^ + xy - 2x ^ 2y + 2) + C^y^ = 0. (29) 

The algebraic curves (1^51) and (P^)) can be parametrized as follows: 

The parameter t — {x + y — 2)/2C — y/2 is undefined if C = and x + y = 2, which is a 
special case of (|29)) . We keep in mind the special case x + y — 2. 

Now we have to check when the elimination of d^F/dx^, d^F/dy"^ from (fTU|) using 
equations pT|) ~ p^ and parametrization ([50)1 leaves the quotient of the coefficients to 
dF/dx, dF/dy equal to y/x, so that we could get rid of the dF/dx, dF/dy terms by 
adding a rational multiple of dF/dt as in ([9]). If a; + j/ = 2, we are led to compare 
coefficients to first order partial derivatives in (|22|) . There remains the case (pO)) with 
C 7^ 0. We replace C i-^ s^, i i-^- t/s and get the parametrization of x and y as in (pO)) . 
The derivative d^F/dx^ in pH)) is reduced by pT|) " P^ to the following expression: 

2cit ^ 4fe2(t + s) 2a+2fei-262-2ci + l\ 9i^ 



f^ -I f^ + 2st+l t J dx 

Abit 2c2{t + s) 2a~2bi + 2b2-2c2 + l\2{t+s)dF f Aabi Aab2 



^2-1 t'^ + 2st + l t + s J s2-l dy yf^-l t'^ + 2st+l 



F. 



To have the right quotient of coefficients to the first order partial derivatives, we must 
have ci = 26i, C2 = 262 and 2a + 1 = 26i + 262. □ 

Now we indicate univariate hypergeometric solutions of the ordinary differential equa- 
tions implied by Theorem [ 



Theorem 2.4 The following function pairs satisfy the same second order ordinary dif- 
ferential equations (with respect to x or t): 



• Fo 



• F. 



a; 61,62 
261,262 



x,2-x] and (a;-2)"''2Fi 



a g+l 
2 ' 2 



61 



(i-t) 



&i + 62 - 5; bi,b2 
26i, 262 

'61 



4t 



{l-s)ist^~l) 



2bi+2b2-l p 
2i^l 



(i- 
f ^^2-^62 

61 + i 



1)2' 

1 

2' ' 



sit~ir 



st^ 



In the last pair, s is any constant. 

Proof. The differential equation ((22|) for the ^2(2;, 2 — x) function has the foUowing 
singularities and local exponents: 



at 


x = 


0, 


1-261, 


at 


x^2 


0, 


1 - 262, 


at 


x = 1 


0, 


61 + 62 - a. 


at 


X — QC 


: a, 


61 +62. 



Local exponent differences at a; = 1 and x = cxd are both equal to 61 + 62 — a. There 
is a chance that (|22p is a pull-back of Euler's equation ([6|) with respect to the covering 
z 1-^ a^^/(2 — xY , with its local exponent differences possibly the following: 

l-C=i-6i, A-B = \-b2, C-A-B = bi+b2-a. 

The shift of the local exponents at a; = cx) and a; = 2 by —a would be needed. Therefore 
we consider Euler's equation ^ with C = l/2 + 6i, A = a/2-|-l/2 — 62, B = a/2, compute 
its puUback with respect to the covering z i~-+ x^ 1(2 — x)'^ and projective normalization 
y{x) 1-^ {x — 2)^"'y{x), and check that the obtained differential equation indeed coincides 
with (f22l) . 

To get the result for the second F2 function, we first consider differential equation 
(1^5)) for the function pp]) of Theorem 12.21 The differential equation has the following 
singularities and local exponents: 



at t = 1 and t = —1 : 

at the roots of t^ + 2st + 1 

at t ~ 00 : 



0, 1-261, 

0, 1 - 262, 

261 + 262, 261 + 262-1. 



If the point t = cx) is an apparent singularity, equation (|23p might be a pull-back of 
Euler's equation ^ with respect to the covering z 1-^ K{t + l)^/(i — 1)^, with the local 
exponent differences 



l-C = i-6i, 



A^B^ 



bi, 



C-A-B = l-2b2 



and the constant K = {s + l)/(s — 1) adjusted so that the other 2 singular points he 
above z = 1. The shift of the local exponents at t = oo and t = 1 by 1 — 2bi — 2^2 would 
be needed. 

Indeed, this is a puUback. We conclude that the function in (|20p and 



^l_ty--2b,-2,,^^^ 







(31) 



satisfy the same second order differential equation. After substituting ti-^(t + l)/(i— 1) 
and s 1-^ (s + l)/(s ^ 1) i^^to the two functions, we obtain the last pair of functions of 
this theorem. □ 

Recall ,B, Section 2.9] that Euler's equation ^ has 24 hypergeometric 2F1 series 
solutions (representing 6 different functions) in general, as discovered by Kummer. Gen- 
eralized hypergeometric equation (jlOOp has 6 hypergeometric 3F2 solutions in general. 
Therefore one can take alternative univariate hypergeometric solutions in Theorems 12. II 
and l2.4l the presented ones look most convenient representatives. In Appendix, we recall 
similar sets of Appell's solutions, say (|102p - (|107p . to the classical systems like pT|) -P ^ . 

Now we discuss explicit relations between the function pairs in Theorems 12.11 and 
12.41 In each case of Theorem 12. 1[ the two functions differ by a constant multiple in a 
neighborhood of x = in general, because they satisfy the same Fuchsian differential 
equation, start with the same local power exponent at a: = 0, there are no other integer 
local exponents at a; = in general, and both functions are defined by proper power 
series around x = 0. In the first case, the two functions are equal to each other, trivially. 
For the second function pair, it is tempting to expand the F2{x,l) function using the 
power series definition ^, and use Gauss' formula Jj Theorem 2.2.2] to evaluate the 
2Fi(l) coefficients to the powers of x. That would give the following formula: 



Ci,C2 



^ I u^„ r(c2)r(c2 -a-62) p / a, 61, a - C2 + 1 

r(c2 - a)r(c2 - 62) \ Ci, a + 62 - C2 + 1 



(32) 



Unfortunately, this identification is incorrect, because the 2Fi(l) coefficients to high 
enough powers of x will eventually divergqj; see [TS] . 

In the first case of Theorem l2.4[ the F2{x, 2 — x) series does not converge. Linear rela- 
tions between F2{x, 2 — x) and 2F1 solutions of the same second order Fuchsian equation 
are cumbersome in general, since analytic continuation has to be identified scrupulously. 

Nevertheless, if the parameters 61 and 62 in the F2{x,2 — x) series are negative inte- 
gers, then the series becomes a finite sum of hypergeometric terms, while 2Fi(x) functions 
have a dihedral monodromy group. The dihedral 2Fi(2;) functions can be expressed el- 
ementarily in terms of nested radical functions, and the terminating F2 series occur in 
the most general form of these elementary expressions. We discuss this in Section [T] 



^In the published version of this paper, in J. Math. Anal. Appl. 355 (2009), pg. 145-163, several 
wrong formulas of this type were claimed. The subsequent paper | 18| rectifies and explains the error. 
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In the second case of Theorem 12.41 let us substitute s i— > s^, i i~^ t/s. We conclude 



that 



and 



Fr 



bi+b2- |; 61,62 
261,262 



4st 



{s'-m'-l) 



(t-s) 



I 61 + i 



2F1 



61 + 62 - 

61 + 



i62 



t' 



(33) 



(34) 



satisiy the same second order ordinary differential equation, with respect to t. In this set- 
ting, the 2Fi(s^) expression is just a constant factor. But the symmetry between s and t 
suggests that equation system pT|) -([T ^ for i^( ''^^''Ih^ib'''^'^^ ^iv) can be transformed 
following p3p to a differential system where the new variables s and t are separated. 
Moreover, the two separated equations with respect to s or t would be simple trans- 
formations of Euler's equation ^. In Section 2] we relate this situation by a known 
transformation to the well-known case ([5]) of separation of variables for the F^^x^y) 
function. 

Here is an identity illustrating the last case of Theorem 12.41 obtainable after substi- 
tution s 1-^ {s + l)/(s — 1) into (P5)) ~ ([M)) and consideration of F2 and 2F1 solutions of 
the same differential system around the point (i, s) = (0, 0): 



61 



l + s + t-st 
1-s 



62 - i; 61,62 
26i,2"62 

2fci+262-l 



4(l-s2)t 



4(1 -t^), 



{1 + s + t- sty ' [i + s + t- sty 



2F1 



61+62-5,62 

26, 



4s 



1 



2F1 



61 + 62-^,6 



61 + 



(35) 



One can check this identity by comparing low degree terms of bivariate Taylor series 
around (t, s) = (0,0). 

3 Identities with Appell's F3 function 

A system partial differential equations for Appell's F3 function is: 



d'^F d'^F dF 

^(1 - ^) a^ + yg^ + (^ - («i + ^1 + 1)^) -d^ - «i^i^ = 0' 

Q2p Q2p gp 

y(^ - y)irT + ^^inr + (^2 - (a2 + 62 -H l)y) ^— - a2b2F = 0. 
oy^ oxoy ay 



(36) 

(37) 



This is an equation system of rank 4. The singular locus of this equation on P"'^ x P^ is 
the union of the following curves: 



a; = 0, X — 1^ X — oo, y = 0, y = 1, y = oo, xy = x + y. 



(38) 



There is a straightforward relation between F2{x,y) and F3{x,y) functions, coming 
from a transformation between differential systems (fTT1) - P^ and ([5Sl) - ([57)) . The relation 
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is that the functions 
a; 61,62 



R 



Cl,C2 



x,y 



1 + 61 + 62 — a 



1 1 

x' y 



(39) 



both satisfy ([TT|) -(fT2 |) . Therefore we can translate the resuhs of Theorems 12.11 through 
2.41 straightforwardly. 

Theorem 3.1 The following function pairs satisfy the same ordinary differential equa- 
tions (or order 2 or 3, with respect to x or t): 



• F^ 



F. 



F3 

„l-c 



fll 


a2; 61, 
c 


b2 


ai 


a2; 61 
c 


b2 


ai 


a2; 61 
c 


b2 



X, I and 2F1 



X, 1 I and 3F2 



and 



ai, 61 



ai,bi,c- 02 - 62 
c— a2,c — 62 



• Fs 



x'-' [1 ~ xr^ 3F2 

1-61,1 -62; 61,62 



X- 1 

1 + fll + a2 — c, 1 + 61 + 02 — c, 1 — 62 
1 + fll + 02 + 61 — c, 1 + 02 — 62 



l~x 



2a; -1 



(1 - xy-' (1 - 2x)'- 2F1 



i(62-6i+c),i(6i + 62+c-l) 



4a;(l-x) ; 



• F. 



1 - 61,1 - 62; 61,62 

3/2 



{t-iy s{t^lf 



At ' (i-s)(st2-i; 



and 



t'^{i-tr'2M\'lf 



st^ 



St2-1 



In the last pair, s is any constant. The function pairs represent all cases when a uni- 
variate specialization of Appell's F3 function satisfies a second order ordinary Fuchsian 
equation fully implied by ([5S)) ~ (|37p . 

Proof. Apart from the first trivial case, we employ correspondence (|39p in each case of 
Theorems \TJ\ and ?U\ We substitute 

xi-^—, yi-^-, ci 1-^ 1 + 61 - ai, 021-^1 + 62 — 02, a 1-^ 1 + 61 + 62 — c, 
X y 

get rid of the power factors x'''^ y^^ , and choose sometimes a more convenient companion 
from the 24 Kummer's or six 3F2 series of the second function. □ 



In the first and fourth function pairs, a direct identity between the two functions 
holds in a neighborhood of a; = 0. In the second pair, the following formula around a; = 
holds if Re(c — 02 — 62) > 0: 

ai,a2; 61,62 



F. 



X, 1 



r(c)r(c-a2 - 62) /ai,6i,c- 02 - 62 

r(c- a2)r(c- 62) V C-02,C-62 



(40) 
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Different from (|5^ . here convergence of the 2Fi(l) coefficients to powers of x in the 
i<3(a;, 1) expansion actually improves with higher powers of a;; see |18) . In the third 
and last cases, relations between Appell's Fj, and 2F1 solutions of the implied ordinary 
differential equations can be derived following their analytic continuation. 

The F^{x,x/{x — 1)) and Fj,{x,x/{2x — 1)) functions are considered in [^. Karlsson 
expresses i^3(a;, x/ {x — 1)) functions with c = oi + 02 + 61 + 62 or c = oi + 62 in terms 
of Gauss hypergeometric functions; these relations can be obtained by simplifying the 
respective 3F2 function in Theorem 13.11 Karlsson's relation between F^{x,x/{2x -— 1)) 
and 2Fi(4a;(l — x)) functions differs from ours by Euler's transformation [2 (2.2.7)]. 

Translating the case ((55)) - (IMl) of variable separation, we conclude that the bivariate 
functions 



and 





\ 


- 61, 1 - 62; 61,62 

3/2 


Ast ' 


(t-s)' \ 




{s^-lW-l)J 


s''if''i(l- 


s-t 


2F1 


2-5,62 


^ 

J 


'^\ b^ + h 



satisfy the same system of partial differential equations with respect to s and t. 

4 Identities with Appell's F4 function 

A system partial differential equations for Appell's ^4 function is: 



a:;(l — x) 



yi^~y) 



d^F 
dx"^ 



d^F 



, d^F 



, d^F 

dx'^ 



2xy 



d^F 



dF 



^ a +{ci-{a + b+l)x) 
oxoy ox 



d^F 



dF 

-{a + b+l)y^—-abF = 0, 

dy 

dF 



dF 

-(a + b+l)x- abF = 0. 

ox 



(41) 



(42) 



(43) 



(44) 



The singular locus on P is x — 0, y — 0, the line at infinity, and 



x^ +y'^ + 1 = 2xy + 2x + 2y. 



(45) 



The quadratic singular locus can be loosely described by the equation ^/x + ^ = 1. It 
can be parametrized as follows: 



a; = t^ 



v^i^-ty 



(46) 



First we characterize univariate specializations outside the singularity curves. The 
general case can be identified as Bailey's case ([8|) of variable separation. 
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Theorem 4.1 Suppose that a univariate specialization of Appell's F^^x^y) function is 
not onto singular locus xy = or GS^, and satisfies a second order ordinary differential 
equation fully implied by (|15)) - (|H)l . Then the specialization is represented by the general 
expression 



Fa 



1 



st, (l-s)(l-<) 



c,a-\-h — c 
where s is a constant. The general second order equation is satisfied by 2F1 



(47) 



a, h 
c 



Proof. The rank of the differential system (|^5)l - (Ml) is 4, and the partial differential 
equations of order 2 are linearly generated by (|15)) - (P^ . Linear elimination oid^F/dxdy, 
d^F/dy^ from expression pI7|) for d^F/dt^ gives the following coefheient to d'^F/dx^: 



■ 2 , -^ y ■ ■ , -^ -2 

X H xy -\ — y 

y y 



(48) 



This expression must be equated to 0. The expression is a quadratic form in the deriva- 
tives, with the discriminant equal to 4:{x'^~2xy+y'^~2x—2y+l)/y'^. We can factorize p8| 
into linear differential forms if we parametrize the surface w'^ = x'^ — 2xy + y'^^2x — 2y+l 
A parametrization is easy if we set w = y + u; then 



y 



(x + u — l)(x — u — 1) 



2{x 



Expression (|48p factorizes as follows: 

(x^ + 2xu + u^ -2x + 2u' 



1)^ 



2{x + u+ l)^{x + u- l){x -u-1) 



{du + dx) [xdu — (m + l)dx) , 



(49) 



(50) 



where we write du, dx instead of du/dt, dx/dt because the variable t is irrelevant. The 
non-differential numerator factor defines the singular locus, since its parametrization 
a; = t^, u = — (f — 1)^ is translated to PS)1 by (|i^. The differential factors give the 



following two families of solutions: 



^Cx-l. 



They translate into the following relations between x and y: 



y 



C-\ 

cTT" 



\-c 



y 



\-c c-\ 



c + i 



(51) 



(52) 



The two solution families actually coincide, since they are related by the parameter 
transformation C 1-^ (3 — C)/(l + C). After the substitutions C i~-> 2s — 1, x i-^ st, the 
first family is described by y — (1 — s)(l — i), giving us (|47p. 
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Adding equation (j43p multiplied by s with equation 
substitutions x = st, y = {s — l){t — I) gives us 



til-t) 



IfC2 






+ 2s(s-l) 



d'^F 
dxdy 



+ (^-iy 



9y2 



multiplied by 1 — s, and the 



+ (ci - {a + h+l)t)s 
dF 



dx 



+ {a + h+l-C2-{a + h+l)t){s-l)—-ahF = 0. 



6 + 1 — ci, we recognize Euler's equation ([6]) in i. 



D 



For fixed s, the arguments x and y in (j47]) are linearly related. The linear relations 
are precisely the tangent lines to the singularity curve (|45p . If we let both s and t vary, 
their symmetry indicates Bailey's case JU of variable separation for differential system 
(gSll-lIlll) when ci+C2=a + 6+l. The conclusion is that 



Fa 



c, a + 6 — c+I 



st, (l-s)(l-i) and 2F1 



a, h 
c 



2F1 



a, 6 



(53) 



satisiy the same system of partial differential equations (with respect to s and t). This 
is illustrated by identity ([5]). 

The F2{x,y) case of variable separation in ([551) - ([M)l is related to Bailey's case via 
the following identity, with c = 2b: 



Fa 



a: b 



x,y 



i^ + yY 



2ap(a-Aa-b+\ 
^\c,2a-2b+l 


X 


4y 


[l + yf 


[i + vY 



(54) 



This identity is derived by Srivastava in [15]; it is presented in [T6l 9.4.(215)]. An equiv- 
alent identity is obtained in [3', pg. 27], and presented in T^, 9.4.(97)]. The Appell 
functions in p3p and (j47p are generic F2 and F^ functions satisfying the quadratic con- 
dition of Sasaki- Yoshida [m Sections 4.3, 5.5]: four linearly independent solutions of 
their partial differential equations systems are quadratically related. 

In [m Section 5.4], the following identity between F^ and F2 functions is implied: 



2' 2 



61- 



2 2 



(l+x + y)-'^F2 



a; 61, 62 
26i,2&2 



2x 



22/ 



y +1 X + y + 1 



(55) 



Equivalent identities are presented in [5], [ISl 9.4.(175)] and [TT, Lemma 5.2]. This 
identity has no application for i^2(a;, 2 — x) functions. 

There remains to consider specialization of the F^ function to the quadratic singularity 
locus (H5t . Due to parametrization (|^5)) . we are looking at Fi{t^ ^ (t — 1)^) functions. Can 
such a function satisfy a second order Fuchsian equation? 



Theorem 4.2 The univariate functions 
a- b 



Fi 



c, a + & — c+l 



i,{i-ty 



and 



2F1 



2a, 2b 
2c -I 



(56) 



satisfy the same Fuchsian equation. This is the only case when the differential system 
\4-S^-^^^ fully implies a second order Fuchsian equation. 
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Proof. Let J be the differential ideal generated by (|33])-([331). After substitution ([i^ 
into the coefficients, the rank of partial derivatives becomes 3, and there is another 
independent linear relation between the partial derivatives of order at most 2. Yet, 
second order partial derivatives are eliminated from the full differential 

already by equations ((431) -(|44]). After we substitute (j46]) into the coefficients in (|43|) - (|44l) . 
add the first equation multiplied by 4/(1 — i) and the second equation multiplied by 4/i, 
we get 

,2d^F o/ -.d^F ,, .2d^F (a + b+l)t-cidF 
4t^^+8i(.-l)^+4(t-lf^+4^-^-^ ^ 

^,i^ + b+l)t + c,-a-b-ldF 4ab (58) 

t dy t{t-l) ^ ' 

Identifying ([57]) and 



dF dF , ,dF 

-.2t- + 2(^-1)- (59) 

we recognize in (jSSp the expression 

d^F {2a+2h+l)t + l-2ci dF 4ci+4c2-4a-46-6 5F 4a6 

"di2"^ t(t-l) It^ t 9^"^ t(t-l) ~ ■ ^ ' 

We have differentiation here with respect to t only if ci +C2 = a + 6+ |. Eulcr's equation 
is recognizable. □ 



t^, (t — 1)^ ) is computed by 



The general ordinary differential equation for i^ 

VCi,C2 

Kato Ji), (4.1)]. We reproduce the differential equation in Theorem 16.11 here below as 
demonstration of general computational techniques; see formula (1901) i'^ Theorem 16.11 
The most important knowledge at this point is that the differential equation is Fuchsian 
with singularities ati = 0, t = l,f = oo, and the following local exponents: 



(61) 



If some difference of local exponents at the same singular point is equal to 1 , then there is 
a chance that the accessory parameter is right and the Fuchsian equation is projectively 
equivalent to Appendix equation pOOp for the 3F2 function. We present these cases in 
the following theorem. 

Incidentally, we recall some Kato's results on reducibility of the monodromy groups of 
differential system (|15)) - (H^ or the ordinary equation for Fi{t^ ^ (l—t)'^). The monodromy 
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at i = : 


0, 


2-2ci, 


C2 — a — 5; 


at i = 1 : 


0, 


2 - 2c2, 


ci - a - 6; 


at t = 00 : 


2a, 


26, 


Ci + C2 - 1 



group of differential system (I43p - (P^ with ci,C2 ^ Z is reducible if and only if one of 
the numbers 



a, b, ci — a, ci — b, C2 — a, C2 — b, ci + C2 — a, ci + C2 — b 



(62) 



is an integer [TUl Section 8] . If no difference of the local exponents in (|CT|) at a singular 
point is an integer, then Kato's ordinary equation is reducible if and only if one of the 
numbers in ([62| or ci + C2 — a — 6 — i is an integer [lOl Section 14] . 



Theorem 4.3 The univariate functions 
a; b 

■-^ 2' 2 



Fi 



t'Ai-tf 



3^2 



2a, 26, c 

+ 6+i2c 



satisfy the same ordinary Fuchsian equation. That equation is satisfied by 

■a+i;6+i 



{l-t)F,\ 



(l-t)-2«F4 



2' '^ ' 2 

1 3 

2' 2 



i^(l-o^ 



a; a + i 



(i-l)2^(i-l)^ 



(63) 

(64) 
(65) 



as well. This represents all cases with c\,C2 ^ Z when the monodromy of (j43P " (|44p is ir- 
reducible and Kato 's differential equation (IffOp is projectively equivalent to hypergeometric 
equation ^10U\i . 

Proof. Here are all cases when a difference of local exponents is equal to 1: 

• C2 = ^. Once we identify the local exponent sets (pTjl and (|10ip under this condi- 
tion, Kato's equation coincides with pOOp up to variable notation and a polynomial 
factor. We get the relation between the functions in ([521) after the substitution 
ci 1-^ ci + i (for a better form). Functions (|64p - ([65)) are solutions of the same 
ordinary equation by Appendix transformations (jl08D - (|110p . 

• C2 = |. To identify the local exponents in fCT]) and poip . we apply the projective 
transformation y{t) i-^ {1 — t)~^y{t). This shifts the local exponents at i = 1 by 
1, and the local exponents at t = oo by —1. We get the coincidence of Kato's and 
hypergeometric equations unconditionally. A direct conclusion is that the functions 



F4 



a; b 

^ 3 



t'Ai-ty- 



and 



1 



1-t 



2a- 1, 2fe- 1, c- 
a + &- i,2c- 1 



t 



satisfy the same ordinary Fuchsian equation. After multiplication of both function 
by 1 — t and a simple transformation of parameters we get the 3F2 function as in 
(j63p . and the Fi function as in ([64]) . 



• ci ~ a + b + 1. The local exponents in (|6ip and (|10ip can be identified directly. 
The accessory parameter is right if C2 = 5 '^'" "^ ~ or 6 = 0. In the later two cases, 
the monodromy group of the system (|15)) - (P^ is reducible. 
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• ci = a + b — I. The local exponents in (pT]) and (jlOip can be identified after the 
projective transformation y{t) i— > (1 — t)~^y{t). The accessory parameter is right if 
C2 = I or a = 1 or 6 = 1. The later two cases are reducible. 

• ci + 2c2 = a + b+ 1. The local exponents in ([CT|) and (jlOip can be identified after 
the projective transformation y{t) h^ (1 — t)^~'^'^^y(t). The accessory parameter is 
right if C2 = -^ or C2 = a or C2 = &. The later two cases are reducible. 

• ci + 2c2 = a + 6 + 3. The local exponents in (|CT|) and (jlOip can be identified after 
the projective transformation y{t) i-^ (1 — t)^^'^'^^y{t). The accessory parameter is 
right if C2 = I or C2 = a + 1 or C2 = fo + 1. The later two cases are reducible. 

• 6 = a + i. To identify the local exponents in ((6T|) and (|10ip . we permute the 
singular points t — \ and t = oo, and shift the local exponents there by ±2a. Then 
Kato's and hypergeometric equations coincide unconditionally. A direct conclusion 
is that the functions 



a, a , 2 

Cl,C2 



' tMl-tf) and {l^tr.YJ'\'^-'^^^'/--'^ 
I \ 2a — C2 + 1,/Ci — 1 



t 



t-1 



satisfy the same ordinary Fuchsian equation. After straightforward transformations 
we get the 3F2 function as in (|55)l . and the F^^ function as in (pS)) . 



• ci + C2 — 2a. To identify the local exponents in (|61|) and (llOip . we permute the 
singular points t = 1 and t = 00, and shift the local exponents there by ±(2a — 1). 
The accessory parameter is right if a = 6 + ^ or a = 1 or ci = a. The later two 
cases are reducible. 

• ci + C2 = 2a + 2. To identify the local exponents in ([CT]) and (jlOip . we permute 
the singular points t = 1 and t = 00, and shift the local exponents there by ±2a. 
The accessory parameter is right if6 = a+^ ora = Oorci =a + l. The later two 
cases are reducible. 

• The cases a = ^, ci = |, C2 = a + 6+1, C2 == a + 6 — 1, 2ci +C2 = a + &+l, 
2ci + C2 = a + b + 3, a — b+ ^, ci + C2 = 26, ci + C2 = 26+2 are symmetric to the 
above. 

Non-reducible cases fall into the cases when ci or C2 is equal to i or |, or a — 6 = ±^. D 

Similarly, we may find the cases when Kato's differential equation is a symmetric 
tensor square of Eulcr's hypergeometric equation ^. Recall that if yi, 2/2 is a basis of 
second order linear ordinary differential equation, then yi, j/|, 2/1^2 form a basis of the 
symmetric square differential equation. One case when Kato's differential equation is a 
symmetric square of ([6]) follows from Bailey's case of variable separation. If we set s — t 
in (1531). we conclude that the functions 



^^1 ^6' 

c, a + 6 — c - 



t\{l-ty] and 2F1, 
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a, b 



t (66) 



satisfy the same ordinary Fuchsian equation. AppcU himself derived this fact exphcitly 
[2j. Here is characterization of other such cases. 

Theorem 4.4 The univariate functions 



Fa 



2c- 



3c -1 



c+\,c- 



t\{l-tf 



and 



2F1 



c, 3c — 1 
2c 



(67) 



satisfy the same ordinary Fuchsian equation of order 3. That equation is satisfied by 

c; 2c 



t 



l-2c 



1 
2 

C, C+ T 



t'A^-ty 



il-2'=(l-t)l-2c^^ 



tMl-tf 



(68) 



(69) 



as well. Together with (j^ffl) this represents all cases when Kato's differential equation 
(|ffOI) is projectively equivalent to a symmetric tensor square of Euler's hypergeometric 
equation ©. 

Proof. A necessary condition on local exponents (pT|) is that at each singular point one 
local exponent is the arithmetic mean of the other two. If we pick for the arithmetic 
mean at i = the exponent C2 — a — 5, or at i = 1 the exponent ci — a — 6, or at i = cxd the 
exponent ci + C2 — 1, we arrive at Appell's case (j66p . Up to symmetries of parameters, 
there are three different ways to pick up other local exponents as arithmetic means: 

• The exponents 2 — 2ci, 2 — 2c2, 2a are the arithmetic means at t = 0, t = I, 
t — 00, respectively. We have 3 linear equations for the 4 parameters, and the 
general solution is presented in (j67p . The local exponents of the 2F1 function must 
be equal to the non- arithmetic mean exponents divided by 2. Explicit computation 
confirms coincidence of ordinary differential equations; in other words, the accessory 
parameter turns out to be right. Functions (|68|) - ((69| are solutions of the same 
ordinary equation by Appendix transformations (jl08D - (|110p . 

• The exponents 2 — 2ci, 0, 2a are the arithmetic means. A general parametrization 
of the parameters is the same as in (|B5| . Evidently, we get a related case. 



• The exponents 0, 0, 2a are the arithmetic means. A general parametrization of the 
parameters is the same as in (j69p . □ 



5 Identities with Appell's Fi function 

A system partial differential equations for Appell's Fi function is: 

Q2p Q2p Qp Qp 

^(1 ~ ^)-^ + y(l ~ ""^d^ + {c-{a + bi + l)x) — - biy— ~ abiF = 0, (70) 

g2p Q2p QP QP 

2^(1 - y) a^ + ^(1 - y^O^j + (c - (a + ^2 + l)y) -g- " b2X— - ab2F ^ 0. (71) 
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This is an equation system of rank 3 generally. The equation 



dF \ { BF 

(1 - y)— - hj [Eq. m] + [i^ - 1)^ + h ] [Eq. dZID] 



can be computed to be a — c + 1 times 

d^F 



(y-x) 



dxdy 



dF dF 

ox oy 



(72) 



(73) 



This equation is of order 2 as well; it holds even for the Fi functions with c — a + \. In 
the case c = a + 1, the equation system ([7D)) - ([7T|) is not holonomic; equation ([75]) has to 
be added to get a finite dimensional space of solutions. 

The singular locus of this equation on P^ x P^ is the union of the following lines: 



a; = 0, X = 1, X = oo, y = 0, y = 1, y = oo, y = x. 
Univariate specializations to the singular lines are known since [3]: 



(74) 



Fi 
Fi 



a; &i,62 

c 
a; &i,62 

c 
a; &i,62 

c 



\ 

X, 
/ 


= 2^*1 


^)^ 








(75) 


x,l\ 


r{c)T{c^a-b2) 
T{c-a)Tic-b2f ' 


a, 6i 
.c-62 


^y 


(76) 


X, X J 


— 2^*1 


+ &2 


^)^ 






(77) 



The second identity holds if Re (c — a — 62) > 0. 

Theorem 5.1 Suppose that a univariate specialization Fi(x{t),y{t)) is not onto any of 
singularity curves in (TO|. Then equations ( | 70^ ~( f7I| ) and ^7M fully imply a second order 
differential equation for Fi {x{t), y{t)) if and only if the functions x{t) and y{t) are related 
by the non-linear differential equation 



- - - - a + 1 

X y \x — I J/— 1 

x{x- 1)2/(2/- 1) fx y 



y 



[y - x)xy 



x{x-l) y{y-l) 
y \(bix 622/ 



X yj\x{x-l) y{y-l)J\x-l y~l 



0. (78) 



Proof. The second order derivative d^F/dt^ in pil)) can be reduced by equations 
([70)1 - (|7T|) and (|73p to first order partial derivatives without additional restrictions on 
x{t),y(t). From equation 



[Eq. ^] 



y 



x{l-x) ■ 2/(1-2/) 

we straightforwardly derive 



[Eq. dUl)] 



(xy - xy)'^ 
xy (2/ - x) 



[Eq. dZSl)] 



d^F 



\. _ c-{a + bi + l)x ^^ ^ b_2X^ ^ b2{x-yf ^ 62(2: + 2/ - 1)2/^ \ dF_ 
x{l — x) X y — X 2/(1 ^ y) J dx 
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.. _ c-{a + b2 + l)y .2 ^ hf _ bi{±-yf ^ bijx + y - l)x'^ \ d£ 
y{l-y) y y-x x{l - x) ) dy 

abix^ ab2y^ 



x{l-x) y{l-y)' 

We can recognize the derivative dF/dt on the right-hand side precisely when the quotient 
of the coefBcients to dF/dx and dF/dy is equal to x/y, according to ([9]). This leads to 
equation dTH]). □ 

Without loss of generality, we can set i = 1, i = in equation (|75|) . Then the non- 
linear equation can be written in the form y = gs^y), where 33(2/) is a cubic polynomial in 
y with coefficients rational functions in x, y. This is not a canonical form of a non-linear 
differential equation without movable essential singularities [U Chapter XIV]. Neverthe- 
less, its Liouville invariants [71 (14)-(I5)] are zero, hence the equation is a differentiable 
transformation x = ^{X, Y), y = ^{X, Y) of the simplest equation Y" — 0. However, it 
is not straightforward to find the transformation. 

Using the identity 



Fi 



a; 61, 62 
c 



x,,) = (i-.xr''^F3|^-«';;=^^''^ 



x~ 1 



,y 



(79) 



one can translate the two non-singular cases of Theorem 13.11 to a few univariate special- 
izations of the Fi function relevant to Theorem 1 5. II The ^2(2:, 2 — x) or i^3(x, x/{2x — 1)) 
case can be translated as follows. 

Theorem 5.2 The following two univariate functions satisfy the same ordinary Fuch- 
sian differential equation of order 2: 



Fi 



a; 26, a — 6 



1 + 6 

The same equation is satisfied by 
'a; 1 — a, a — 6 



and 



(l-x)- 



2^1 



4.T 



{x-iy 



(l-2:)-"i^i 
x-^Fi 



1 + 6 
a; 26, a — 6 
2a 
a; a ^ b,a — b 
1 + a - 26 



x-1' 



1 — X, 1 — x' 



-, X 



x--F^ ' 
{x-\)-''Fr 



x-\\ - xy-^^Fi 



a; 26, 1 — a 

1 + 6 

a; a — b,a — b 

2a 
a; 1 — a, a — 6 
^ ' 1 + a - 26 
1 — 26; a — 6, a — 6 
1 + a - 26 



x-'il-xY 



'F 



1 — 2b; a — 6, 1 — a 
1 + a - 26 



X 


x^ \ 


X+V x^-lj 


.T-1 \ 


,1-a; 


X J 


^ XI 1] 




l-x J 


1 X \ 


l-x'x- ij 


— a 


X + I, X 

) 



(80) 

(81) 
(82) 
(83) 
(84) 
(85) 



21 



Proof. The F3 function in ([7^ can be directly related to the penultimate case of 
Theorem 13. II if a = 1 — 62 and c = 2 — &i — 62- The direct conclusion is that the functions 



F, 



1 - 62; bi,b2 
2-61-62 



1 



and 



(1 



1 



1-61, 



'f^i2'-i,;':i 



4.T 



(x-l)^ 



(86) 



satisfy the same second order ordinary Fuchsian equation. We rename 61 1— > 1 — a, 
62 i-^- a — 5, and adjust the power factor to the 2F1 function, so to get the second function 
in (I80|) . As presented in Appendix formulas ()llip - (ll22p . there are generally 60 Fi bi- 
variate series giving solutions of the differential system (|70|) -([7T]l. Investigation of those 
60 Fi series under the given specialization gives the Fi functions in (|5D|) - ([55)) . □ 



The Fi functions in ((80)) - ((85)) represent all possible relations between the 2 variables and 
the parameters among the specialized 60 Fi series, obtainable from Appendix formulas 
piip - (|122p and formally giving univariate solutions of the same ordinary differential 
equation. In particular, the -Fi function in (j86p has the same relation between the 
arguments and parameters (up to permutation of the two arguments and the 61 , 62 pa- 
rameters) as the first function in ((8T|) . The four functions in ([80)l - (|8T|) can be identified 
in a neighborhood of a; = 0, while the two functions in (I82p can be identified in a neigh- 
borhood of x ~ 1. Notice that relation between the two arguments in ([M|) or ([55]) is 
linear. 

The variable separation case (|4T|) - (|42p for the F3 function translates into a reducible 
transformed system ((70)) - ([7T|l . A representative conclusion is that the differential system 
for 



Fi 



2' ' 2 

3/2 



i-s-t)^ 



(s + t)2' (s2-l)(i2_l) 



(87) 



has the following elementary solution: 



1 



{s + tf''{i^s'y-''{i~t^y''. 



If differential condition (jTS)) of Theorem 15.11 is tried for linear solutions y{x), the 
following cases are found: 



• Reducible Appell's functions Fi 
a; 0,62 



a: 61, 62 



and Fi 



a+l 
X, s \ , where s is a constant. 



x,sx ] , Fi 



a; 61, c — 61 



1 — x.l — sx 



• The cases (l84l) . (|85|) with y — 1 — x, y = x — 1, respectively, and a symmetric case 
with y — X + 1. 



The trivial identity Fi 



-1; bi,b2 



c — bix 



0. 
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6 Computational aspects 

Identities between bivariate and univariate hypergeonietric series are usually derived by 
methods of series manipulation, or using integral representations. The method of relating 
differential equations is usually considered as tedious and computationally costly. In [161 
pg. 314] Srivastava and Karlsson characterize this method as follows: 

This method (though simple in theory) is rather laborious in practice, and is 
not very useful for discovering new transformations. 

They immediately mention that Appell himself 2 showed that the functions in (|66p 
satisfy the same differential equation of order 3. 

Here we demonstrate the basic computational routine to translate linear differential 
equations in partial derivatives to ordinary differential equations for univariate special- 
ization of their solutions. In our terminology, we seek to recognize partial differential 
forms of ordinary equations within specialized forms of given systems of partial differen- 
tial equations. 

In the specialization setting x i— > x{t), y i— > y{t) of Definition 11.11 let J denote 
the left ideal of the Weyl algebra C[x,y]{d/dx,d/dy) determined by a given system of 
partial differential equations. We work in the C{t){d/dt) left module M that is linearly 
generated over C{t) by the full derivatives 1, d/dt, d^/dt'^, etc., and the partial derivatives 
d/dx,d/dy,d^/dx^,d^/dxdy, etc. We let d/dt act from the left following the Leibniz 
rule, with the action on the partial derivatives defined by the identity 

d . d . d ,„„, 

dt=^0-x^'d-y' (^^) 

compatible with ([9]). Let P be the submodule of M generated by d/dt — x d/dx — y d/dy. 
The elements of P give expressions of the full derivatives d/dt,d'^ /dt'^, etc., in terms of 
partial derivatives, like the series of formulas starting with ©-(HU]). 

Let J denote the submodule of M generated by special forms of the elements in J. Let 
P denote the union of J and P, as a submodule of M. We propose to use Grobner basis 
computations for P with respect to an ordering of full and partial derivatives suitable 
for elimination of the partial derivatives. 

Here we briefly demonstrate this strategy by computing Kato's differential equation 
[TUl (4.1)] for the general i^4(i^, (1 — t)^) function. We used this equation in Section HJ 
and presented its singularities and local exponents in (|6ip . 

Theorem 6.1 The univariate function i^ 
equation 



t ,(t— 1) j satisfies the differential 



d^F 
dt^ ^ 
/(2ci- 


/a + fo + 2ci -C2 + 1 a + b-ci + 2c2 + l\ d'^F 
\ t ' t-1 ) df^ 
l)(a + 6-C2 + l) ^ ^a + 6-ci - C2 + 1 + 2a6-|-2ciC2 


V 


t2 ' " t{t- 1) 
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i2c2-l){a + b-ci + l)\dF , 2abi2{ci + C2-l)t-2ci + l]_^ ^ ^ 



(t-i)2 ; dt f^{t-i) 



2 



Proof. In the course of the proof of Theorem 14.21 we basicaUy derived equation ([60)1 as 
an element of P. We differentiate ([60|l with respect to t, keeping in mind that d/dt acts 
on the partial derivative dF/dy following ([89| . or more specifically, following ((59|l . We 
get the partial derivatives d^F/dxdy and d^F/dy^ in the computation; we can eliminate 
them using original equations ((i51) - (|i^ . but then first order derivatives dF/dx, dF/dy 
occur. For linear elimination of these derivatives, we use (|59|) and the equation ([60]) itself. 
Equation ([90|) is obtained as 



^ (c.-a->,-l)(2o.+2c.-2„-2l,-3) ^ (.. + .. + l),+c.-a-l,-2 

r(i: — 1)^ t[t — 1) 

where the second term is added to eliminate second order partial derivatives, the next 
term — to eliminate dF/dx^ and the last term — to eliminate dF/dy. □ 

Kato himself [TU] derived equation ([50)1 via a conversion to Pfaffian systems. The spe- 
cialization transformation of Pfaffian systems is particularly elegant, but Kato's method 
requires conversion of a matrix differential equation to an ordinary differential equa- 
tion. The specialization problem for partial differential equations can be considered as 
a restriction problem for ZJ-modules |13j . However, available implementations of algo- 
rithms for D-modules worked very slowly on our examples. Border bases [H] rather than 
Grobner bases look promisingly suitable for our strategy. Towards the end of Section 1 we 
mentioned the simplest computational view of looking for a linear relation between full 
derivatives (in the specialized variable) expressed in terms of a basis of partial derivatives. 

7 Application: Dihedral hypergeometric functions 



a, 



In the first case of Theorem 12 .41 we have a general function of the form 2F1 

The corresponding differential equation is a general pull-back transformation of Euler's 
equation ^ with respect to a degree 2 covering ramified above two (of its three) sin- 
gular points. The pull-backed equation has 4 singularities in general; it has 3 (or less) 
singularities when a local exponent difference below a ramification point is equal to i 
Theorem 12.41 says the general pull-back equation has solutions expressible as an F2{x,y) 
function with x + y — 2. As we mentioned, the F2 series does not converge unless it 
terminates. 

An interesting case is when the parameters 61, 62 in the first case of Thcorcm l2.4l are 
zero or negative integers, say 61 — —k, 62 = —i- The immediate conclusion is that the 
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functions 



2F1 



a g+l I // 
2 ' 2 ~^ 



and (l + Vi)-"i^2 



-2k,- 



2^ 2 

1 + yi'i + Vi 



(91) 



and satisfy the same second order ordinary equation, with respect to z. The monodromy 
group of the 2Fi(2) function is a dihedral group; the local exponent differences at z = 
and z = 00 are the half- integers fc + ^, ^ + \, respectively. Functions with a dihedral 
monodromy group can be expressed as elementary functions, since after a quadratic pull- 
back transformation (ramified above z = 0, z = cx), in our case) the monodromy group of 
the pull-backed equation is a cyclic group, and hypergeometric solutions become simple 
power or logarithmic functions. An example of such an elementary expression is 



2^:^! 




(i-Vi)"" + (i + V^)"" 



Relatedly, 




(1 - Vi)-° - (1 + Vi)-° 

2a ^/z 

log(l -H ^/i) - log(l - yi) 
2^z 



(« 7^ 0), 



l + x/l~7 



(92) 

(93) 
(94) 
(95) 



The 2F1 function in (j9ip is a contiguous family of Gauss hypergeometric equations with 
the same monodromy group. 

The relation between the two functions in ([9T|) allows us to generalize formula ((92)) . 
The generalization is 




(1 + V^)- 


a 


2 


a 



F, 



F, 



a; — fc, —£ 
-2k, -21 

a; — fc, —£ 
-2k, -2i 



2y/^ 2 

1 + ^/i'TT 

271 2 



(96) 



The F2 series on the right-hand side are finite sums with (fc -I- 1){£ + 1) terms. If 
(^ii) 7^ 0, this is an expression of the 2F1 function as a linear combination of two 
explicit solutions of the same differential equation. The identity can be proved by using 
the symmetry with respect to the conjugation of y/z and checking the value of both 

' a,-e 

-2£ 
algorithm. See [17] for details. 



sides at z = 0, which leads to evaluation of 2F1 



2 ) obtainable by Zeilberger's 
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Similarly, a generalization of ([M)) is 



(^: 



k\2)k+i+i ( i\fc^fe+i_-p^( ^- + A;, I + /c + ^+ 1 



12/*: \2)k+l V2J 



(~l)^Z^+^2Fl| 



f+fc 



2 ^V-2fc,-2^ 



2Vi 2 

2Vi 2 



(97) 



To show this identity, one can verify that both sides satisfy the same recurrence relations 
with respect to k and £, and check the identity for a few values (fc,£). 

Recalling the relation between F^{x,x/{2x — 1)) and 2F1 functions in Theorem 13.11 
Karlsson's identity jTB), 9.4.(90)] gives the following generalization of ((95)) : 



I'Pi 



a-i a+£+l 
2 ' 2 



1 + VT^ 



— a — k 



(l-z)^/2x 



k + i,e + i;-k,-e 

a + k + l 



2^/Y^~z ' 2 



(98) 



Note that the Fj, sum is finite for any integers /c,-^, since it is invariant under the substi- 
tutions k H^ — fc — 1 and £ 1-^ — £ — 1. The F2 sums in (|96| and ([97]) can be written as F-^ 
sums, or vice versa in (j98p . by reversing the double summation in both directions: 



Fi 



x,y 



(2fc)!(2Z)! "^^ ^^^ l-a-k-l 



~2k,-2i 
These explicit expressions are analyzed more thoroughly in |17| 



1 1 

x' y 



(99) 



8 Appendix 

Here we recall a couple of relevant facts about ordinary Fuchsian equations, and review 
related explicit results on Appell' functions existing in literature. 

can be written 



The third order hypergeometric differential equation for 3F2 1 



x^{l-x) 



(Py{x 
dx^ 



x{D + S + 1 - (A + B + C + 3)a 



A,B,C 
D,E 

^(Pyix)_ 
dx'^ 



{DE- {AB+AC+BC+A+B + C+\)x 
The singularities and local exponents are: 



dy{x) 
dx 



ABCy{x) == 0. (100) 



at x = : 0, 1 - D, l-E\ 

atx^l: 0, 1, D + E-A-B-C\ 

aX, x — Qo : A, B, C. 



(101) 
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Recall that a third order Fuchsian equation with 3 singularities is determined by the local 
exponents and one accessory parameter. Say, one can add a scalar multiple ofy{x)/{x—l) 
to the left-hand side of equation (jlOOp without changing the local exponents in (|10ip . 

Solutions of differential systems ((n|)-(ini), ^-§7^, l^~^, ^-^ in terms 
of bivariate hypergeometric series were considered by many authors [6, pg. 222-242], 
[Hi pg. 291-305], starting from Le Vavasseur, Appell, Kampe de Feriet, Borngasser, 
Erdelyi, Olsson. 

In general, each Appell's F2{x, y) function can be represented by four F2 series at the 
origin (0,0): 



a; 61,62 

Cl,C2 



x,V 



= il-x)--F2 
= il-x-y)--F2 



a; 


c\ - 


-h 


h2 




C\ 


, C2 




a; 


&i, 


"2 — 


b2 




C\ 


C2 





X 



y 



x — 1 ' I — X 
X y 



1 



y y 



1 



a;ci-6i,C2-62 

Cl,C2 



x+y—\ x+y—\ 



(102) 
(103) 
(104) 



Note that the parameters of the ^^2(2;, 2 — x) function in (PO)) are invariant under these 
transformations; this invariance is discussed in [Ml Section 10]. Besides, in general there 
are four distinct F2{x. y) functions that are local solutions of ([TT|l - ([T2[) at the origin. The 
other 3 functions are represented by the series 



^^-'''v^-''^F2 



x^-^'F2 


^1 + a-ci; l + 5i -ci,62 


y 2-Ci,C2 


y'-''F2 


^l + a-c2; 6i,l + 62-c2 

V Ci,2-C2 


a~ci- 


C2; 1 + 61 — Ci, 1 + 62 — C2 




2-ci,2-C2 



x,y 



x,y 



x,y 



(105) 
(106) 
(107) 



The F2 system pT|) - (IT^ has Horn's H2 series local solutions at {x,y) = (0,oo) and 
{x,y) = (00,0), and four distinct F3 local solutions at {x,y) — (oo,cxd). The latter are 
obtainable by applying relation (|55)) to the four series in (|102p - (|104p . Application of 
relation ([M]) to (|105P " (|107p does not give different F3 series, but only shows invariance 
of ^ under the permutations oi <-> 61 and a2 <-*■ 62 of upper parameters. Hypergeometric 
solutions of the F3 system (P5)) - (|57)) are described similarly. 

The system (|43)) -(|44 | has four generally different F4 local solutions at {x,y) — (0,0). 
Besides ^ we have 



c^-^'Fi 



a — ci; 1 + 6 - 

2- Ci,C2 



Ci 



x,y 



A — Ci ^,1 — C2 



Fa 



V Ci , 2 - C2 


^C2 


2 + a -ci -C2;2 + 6-ci ^ 


- C2 


2-ci,2-C2 





x,y 



x,y 



(108) 
(109) 



Similar sets of four local solutions exist at {x,y) — (0,oo) and (00, 0). A connection 
formula between ^ and the following two functions is presented in [3, pg 26] and [161 
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9.4.(69)], for example: 



y^'^F^ 



a; 1 + a - C2 
ci, 1 + a — 5 



X 1 

y' y 



ci, 1 + — a 



X 1 



(110) 



The full set of Fi solutions to the system (fTOjI -lfTl ]) consists of 60 Fi series in general. 
First of all, the Fi solutions are identified in sextets: 



Fi 



a; 61,62 



x,y 



il-x)-'^il-y)-'-F, 



c- a; 61, 62 



X y 



X — 1 y — 1 



(1 - x)--Fi 



(i-y)-"^i 



a; c- 61 - 62,62 



X X ~ y 



a; 61, c- 61 - 62 



{i-xy-'^-'^ii-yy-Fi 



X — 1 X — 1 
x-y y 
'^-y' y-^ 

c— a;c-6i — 62,62 



[i-xy^^i-yr 



— a — b'- 



F 



c—a; 61, c— 61 — 62 



, (111) 
(112) 
(113) 



=^,f^),(114) 
1-2/ 

x-y 

x-l 



,y -(115) 



Consequently, there are 10 different Fi solutions of ([7n|) - ([7T|) in general. The other 9 are 
represented by the following series: 



Fi 



a; 61, 62 
l + a+61 + 62 — c 

a; 1 + a — c, 62 



i-.,i-y), x-'^^y-^^F,[^-;\Ti~/^^^'^'' 

1 + 61 + 62 - a 



x-^Fi 



1 + a - 61 



i y. 

x' X 



1 + 61 + 62 - c;6i,6; 
61 + 62 - a 
a; 61, 1 + a — c 



y-^R 



1 + a - 62 



1 1 

x' y 
X 1 

y' y 



(116) 
(117) 



{l-x)-''^{l-yY-''~^^Fi 



i-y 

l~x'' 



{l-xr-'^-'^il-yy-F^ 



62, 



^.-blybl-c+l 



Fi 



c — a; 61, c— 61 — 62 

c — a — 62 + 1 
c- a;c- 61 - 62,62 

c — a — 61 + 1 
1 + 61 + 62 — c; 61 , 1 + a — c 
2 + 61 - c 



,1-y , (118) 



1 — a:, 



l-x 



(119) 



-,y , (120) 



6,_e+i -62 p / 1 + 61 + 62 - c; 1 + a - c, 62 



„bi+b2-c( 



ix-yy-'''-'''ii~xy-''-'Fi 



2 + 62 - c 



1 — 61; 1 + a—c, c— 61 — 62 
2-61-62 



X,- 

V 



x — y x — y 
x—1 ' X 



(121) 
(122) 



Besides, Appel's Fi functions can be expressed as F^ series following formula ([7^ . Ap- 
pell's F2 series realize Fi functions as well; the following identity hold in the neighborhood 
of the point {x/y,y) = (1,0) of the blow-up at {x,y) ~ (0,0): 



F 



a; 61, 62 



.,.)=(./.)^^i^.(V!f"" 

' 61 + 62, c 



x~y 



X 



,y 



(123) 
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In the remainder of the Appendix, we summarize identities in [TBI Section 9.4] that 
relate different Appell's functions to each other or to Gauss hypergeometric functions. 
The expressions for Appell's functions with reducible (up to possibly a quadratic trans- 
formation) monodromy groups are formulas [H 9.4, (87), (88), (98), (99), (108)-(114)]; 
several of these formulas are attributed to Bailey or to fS[. It looks convenient to sub- 
stitute X ^ x/{x - 1), y 1-^ y/{y - 1) in Bailey's formulas [M 9.4, (110)-(113)]. We 
mentioned special relations ([54)) - (|55|) between F2 and F4 functions; another similar rela- 
tion is [m 9.4, (216)]. The following 5 formulas are [M 9.4, (190)-(192), (149), (179)]; 
they give a taste of what can be expected from an exhaustive research of univariate 
specializations of Appell's functions to 4F3 or 3F2 functions. The first three formulas are 
attributed to Bailey, and (|127p is due to Burchnall: 



F, 



F. 



a; b, b 
c, c 

a; 61, 62 
2&i,262 

a, a; b, b 



Fa 



Fi 



a; b 

Cl,C2 

a; b 






= 4F3 



4F3 



= 4F3 



4i:'3 



— 4F3 



a a+1 5 ^ 5 
2 ' 2 ' ' 

' 2 ' 2 


^ 




a a+1 61+62 61+62 + 1 
2 ' 2 ' 2 ' 2 


bi + ^,b2 + \,bi +62 


„ ;, a+b a+b+1 
' ' 2 ' 2 


xA, 
J 


Ci,C2,Ci -f C2 - 1 


4a; 


a a+1 b 6+1 
2' 2 '2' 2 

„ c c+1 

•-'2' 2 


-4 


\ 

x' 
J 



(124) 



(125) 



(126) 



(127) 



(128) 
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